Integrable lattices and their sub-lattices: 
from the discrete Moutard (discrete Cauchy-Riemann) 
4-point equation to the self- adjoint 5-point scheme. 

A. Doliwa*, P. GrinevicbJ, M. Nieszporskif and P.M. Santini^ 



Abstract 

We introduce the sub-lattice approach, a procedure to generate, from a given in- 
tegrable lattice, a sub-lattice which inherits its integrability features. We consider, as 
illustrative example of this approach, the discrete Moutard 4-point equation and its 
sub-lattice, the self-adjoint 5-point scheme on the star of the square lattice, which are 
relevant in the theory of the integrable Discrete Geometries and in the theory of dis- 
crete holomorphic and harmonic functions (in this last context, the discrete Moutard 
equation is called discrete Cauchy-Riemann equation). We use the sub-lattice point of 
view to derive, from the Darboux transformations and superposition formulas of the 
discrete Moutard equation, the Darboux transformations and superposition formulas of 
the self-adjoint 5-point scheme. We also construct, from algebro-geometric solutions of 
the discrete Moutard equation, algebro-geometric solutions of the self-adjoint 5-point 
scheme. We finally use these solutions to construct explicit examples of discrete holo- 
morphic and harmonic functions, as well as examples of quadrilateral surfaces in R 3 . 

1 Introduction 

One of the most important methods to generate integrable equations is to start with a 
fairly general integrable system and apply to it systematically symmetry reductions. This 
approach was, for instance, applied to the multicompontent Kadomtsev-Petviashvili hier- 
archy ^5] and to the self-dual Yang-Mills system [33]. On the level of discrete equations, 
it was systematically used, for instance, to generate various integrable reductions of the 
multidimensional quadrilateral lattice equations (see, for example, [S]). 

In this paper we propose a different, but equally relevant, procedure to generate inte- 
grable lattices. This procedure consists in constructing, from a given integrable lattice, a 
sub-lattice which inherits the integrability features of the original lattice. We consider, as 
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illustrative example of this approach, the 4-point difference equation 

4>m+l,n+l — <Pm,n = 5m,n(0m+l,n ~~ <Am,n+l)i (1) 

where g m>n and (j) m>n are functions: 1? — > C, and its sub-lattice, the self-adjoint 5-point 
scheme: 



on the star of the square lattice, where a^ u , c M)l/ and are functions: Z 2 — > C, and 
where the standard notation / mj „ = f(m,ri) is often used throughout the paper. 

We exploit the above fact deriving several properties of the 5-point system © from the 
corresponding (and simpler) properties of the 4-point system 

As it was shown in [27] (see also (7j), equation Q is an integrable discretization (i.e., a 
discretization possessing Darboux transformations (DTs)) of the Moutard equation 

*,uv = F$ (3) 

(the symbol f u denotes partial differentiation: f u = df/du). As it was shown in |24j . 
equation (J2J) is an integrable discretization of the elliptic (if AB > 1) equation 

(A$ u ) tV + (B$ v ) tV = F$. (4) 

The paper is organized as follows. In section [21 we establish the connection between the 
above difference equations (pQ) and (j2j). In sections 01 and we use the above connection 
to derive the Darboux transformation (DT) and its superposition formula of the 5-point 
scheme (J2J (in all its most distinguished gauge equivalent forms) from the the Darboux 
transformation and its superposition formula of the 4-point scheme (0). In section[2]we first 
construct algebro-geometric solutions of the 4-point scheme Q and then the corresponding 
solutions of the 5-point scheme (J2J). Due to the interesting applications of equations (fj| 
and © (discussed in the second part of this introduction), these explicit solutions provide 
examples of discrete holomorphic and harmonic functions and, at the same time, they 
generate quadrilateral surfaces in M 3 . 

Let us devote the remaining part of the introduction to the presentation of the interesting 
applications of equations (JJ) and (j2J) in differential and discrete geometry, and in the discrete 
holomorphic function theory. Equation Q has a very long history. Bianchi constructed the 
superposition formula for the Moutard equation (J3J) in the following form [2] 

$(12) _ $ = ^!($(D _ $(2h (5) 

a 

where $( 1 \ $( 2 ) and <E>( 12 ) are four different general solutions of four different Moutard 
equations (JHJ), while 9 , 9 2 are particular solutions of the Moutard equation satisfied by 
and a is a potential connected to 9 , 9 2 via the conditions 



1% ~ 2 0\, (6) 

1 t> v — t> V v . 



(7) 
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Figure 1: The Bianchi permutability diagram for the Moutard equation. 



Using the key observation contained in ^5], such superposition principle was interpreted 
in [27J as an integrable discretization Q of the Moutard equation: 

$ -» 4>m,n, -» 0m+l,n, $ (2) -> #m,n+l, $ (12) -» &n+l,n+l- 

In addition, DTs and their superposition principle for the discrete Moutard equation (JI)l 
were also constructed in [27) . 

Solutions of the (real) Moutard equation © allow one to construct, via the Lelieuvre 
formulas [2j, asymptotic nets on hyperbolic surfaces from their normal vector. Analogously, 
solutions of the (real) discrete Moutard equation (JTJ) characterize normal vectors of asymp- 
totic lattices Ej, which can be constructed via the discrete analogue of the Lelieuvre 
formulas [5J- 

Also for the elliptic version of equation (fHjl. the well known Schrodinger equation 

*,« + = ™ (8) 

(which is equation Q with = £? = 1), there exists a geometric interpretation within 
the theory of isothermally conjugate nets and the corresponding analogue of the Lelieuvre 
formulas The 5-point scheme (|2|). introduced in |2I] as an integrable discretization of 
equation (jlj), describes the normal vector of a quadrilateral lattice in M 3 , and this embedding 
is obtained via a suitable generalization of the Lelieuvre formulas [23] ■ 

It is remarkable that equations Q and (J2J), which are basic equations in the recently 
developed theory of the Integrable Discrete Geometries, are also the basic equations of an 
integrable discretization of the theory of holomorphic and harmonic functions. Indeed the 
equations 

VVn+l,n+l — VVn,n = H^m+l,7i ~~ VVra,re+l), (9) 

VWi," + + Vv,^+i + ^V^-i = 4^1/! (10) 

undressed versions of equations Q and (J2J), were introduced as basic objects of the discrete 
holomorphic and harmonic function theory in }13(l9|: they correspond to a natural discretiza- 
tion of, respectively, the d and Laplace operators on the square lattice, and the connection 
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between these two schemes was an important ingredient of the theory. In this theory 
was generalized to rhombic lattices; on the level of the 4-point scheme, it corresponds to the 
discrete Cauchy-Riemann equation Q with a nontrivial pure imaginary potential g m ,n- in 
|21| this theory was extended to discrete Riemann surfaces and, in particular, the 4-point 
scheme Q and the affine form 

a M) „(Vv t +i,i> - ipw) + a M _i^(-0 M -i,^ - i)p, v ) + b^viip^v+i - ip^ u ) + ^^-1(^,^-1 - ipw) = 

(11) 

of the self-adjoint 5-point scheme (0), which corresponds to a discretization of the harmonic 
equation, were studied on the Riemann surface. A theory of Dirac operators and the Green's 
function for planar graphs with rhombic faces were constructed in JH]- The connection 
between the above discrete complex function theory and the discrete complex function 
theory [22] based on the integrable cross-ratio equation [2E] was established in [3], where a 
multidimensional generalization of the discrete complex function theory was also considered. 

We remark that there exists another natural integrable discretization of elliptic operators 
on the plane, based on the self-adjoint scheme on the star of a regular triangular lattice |29| . 
resulting in a 7-point scheme. One of the important features of this approach is that this 
operator is factorized in terms of 3-point operators living on the same lattice. A discrete 
complex function theory based on these operators, including the discrete analogue of the 
Cauchy kernel, has been recently developed in 

2 The 5-point scheme as a sub-lattice of the 4-point scheme 

Consider a lattice, i.e. a map x : T> — > V from a grid T> to a linear space V satisfying 
a certain equation E[x] = (the lattice equation). Consider a subgrid T>' C T>\ if one can 
construct, from the original lattice equation E[x] = 0, a new equation £"[03'] = 0, where x' 
is the restriction of x to the sub- grid T>'\ x' = x\x>' T>' — > V, then x' is a sub-lattice of x 
and E'[x'] = is the associated sub-lattice equation. 

Suppose that the original lattice (equation) be integrable; i.e., suppose that one can 
associate with it linear transformations enabling one to construct solutions from solutions 
(Darboux-type transformations), whose superposition exhibits permutability properties (the 
Bianchi permutability diagram). Since the infinite class of solutions generated in this way, 
once restricted, are also solutions of the sub-lattice, this sub-lattice (equation) will also 
be integrable. We remark that the dimensional reductions are simple examples of such a 
construction. 

2.1 Prom the lattice to the sub-lattice 

In this section we start with the general 4-point scheme 

(12) 

where a, j3, 7 and 6 are C-valued functions on T> = Z 2 , and we explore the possibility to 
construct, from it, a sub-lattice defined on the subset T> = H? e (the even grid) of Z 2 consisting 
of points (to, n) such that to + n is even. 
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Proposition 1. The general 4-point scheme Uty) reduces to a 5-point scheme on the even 
grid if and only if the coefficients a, /?, 7 and 5 satisfy the constraint 

Pm,n^hn— l,n&m,n— lTVn— l,n— 1 — 7m,n^m- l,n&m,n— iPm— l,n— 1- (1^) 

Proof. Consider the 4-point scheme (|12j) in the elementary square 

<5m,n = {("i, n), (m + 1, n), (m, n + 1), (m + 1, n + 1)} 

and in its neighbouring squares Q m -i )n -i, Q m -\,n and Qm,n-i, and move to the LHS the 
terms in which 9? is evaluated on the odd grid and to the RHS the terms in which ip is 
evaluated on the even grid: 

fim,n'*Prn+l,n ~\~ ^m^n^m^n-^X — < ^m,n'/'m+l,n+l "771,71 V77i,r?,; 

(%m— 1,71^771,71+1 &m— l.nVm- l,n — Pm—l.n^m.n 7tti— l,nV7Ti— 1,71+1; (14) 

(%m,n— iym+l,n "I - ^m.n-l'/'m^-l — 777i,n— 1^771,71 0m,n— IVtti+1,71— 1; 

7tti— l,n— l^m— l,n Pm—l,n—l^Pm,n—l — ftm— l,n— l^m,?! "m— 1,71—1^771— 1, 71— 1- 

In order to construct a lattice equation not involving <p on the odd grid, it is sufficient 
to impose the zero determinant condition for the system (|14j) . which coincides with the 
constraint (fl3|) . □ 

Remark. The row vector solution v of the homogeneous version of the system (j!4jl reads: 

Vm,n = {'~Yin—l,n—lOtm—l,nCim,n—li 7?n,n7?Tt— l,n— lC^m,n— 1 j . , 

fl A \ ( 15 ) 

C^m— l,nHm,n'Jm— l,n— 1 > C^m,™— l^m— l,n77(i,n J 

and the wanted sub-lattice equation is the solvability condition («J mi „,F mi „) = of the 
Kronecker - Capelli theorem, where 



■F m ,n — {dm,n ^Pm+l,n+l ~i~ ^m^^Pniyni 7tti— l^V'Tn— l,n+l ~i~ Pm—l^n^PmyTii 
ftm,n—l t Pm+l,ri—l ~\~ 7m,n— lVm.n; <^m— l,n— l^Pm— l,n— 1 "I - C^m—l,n~l^Pm,n) 



T (16) 



is the column vector representing the inhomogeneous term of the algebraic system. This 
solvability condition reads: 



l ~1m,ri~1m— 1,71— l^m,n~ 1 (7m— l,n Vm— 1,7J+1 "I - Pm— l,7i^77»,>i) (17) 
^m- l,n/5m,Ti777i— l,n— 1 \0m,n— 1 ^Pm+l,n— 1 "I - r lm,ri—l ( ftm,n)~\~ 
^m,n— l^m— l,n'Ym,n(<^m— l,n— l^m— 1,71— 1 ~\~ Ctm—X,n—'i. t Pm,n) — 0' 

Proposition 2. T/ie /owr poini scheme M/Jjl with non-vanishing coefficients a, (3, 7, <5 satis- 
fying constraint hl'J\) is gauge equivalent to the discrete Moutard (discrete Cauchy-Riemann) 
four point scheme £7J) whose corresponding five point scheme is the affine self-adjoint five 
point scheme Ml)) . 

Proof. Let the potential p : 1? — » C be a solution of the system 

^771,71^771+1,71+1 — Om,nPrn,m ,nPm+l,n — 1m,nPm,n+li 

(18) 
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whose compatibility condition is the constraint (jL-ij) . Then the function ip : Z — ► C 

/Pm,n, (19) 



satisfies equation 
with the potential 



^m+l,n+l — ^m,n — ifm,n(lpm+l,n ~ ^m,n+l), (20) 



_ ■ f _ Pm,nPm+l,n _ '~fm,nPm,n+l \ 
9m,n — 1jm,n — 7 — 7 j V / 

Om,nPm,n Om,nPm,n 



and equation (fTT)) reduces to the 5-point scheme 

Q"m,n(' l Pm+l,n—l 1pm,n) ®m— 1,71+1 \lp7n— l,n+l VVr^n) - ! - 
bra,n{^m+l,n+l ~ ^m.n) + &m-l,n-l [1pm-l,n-l ~ 1prn,n) = 0, 

with the functions a m)n and 6 mjn defined as follows 



(22) 



Q-m,n — fm,n— li ^m,n — • (23) 

Jm,n 

The connection between (j'2U|) and (|22[) we have established here can be easily verified directly 
considering the 4-point scheme ()2(J|) for the vector tp £ V in the elementary square Q m ,n 
and in its neighbouring squares Q m _i ;n _i, Q m -i,n, and Q m ,n-\'- 



(24) 



^-^(V'm+l.n+l — V>m,n) — i(lpm+l,n ~ VVi.n+l), 

fm-i n—i C V^rri — 1 ,ti — 1 ^m,n) — ^(^m— l,7i VVn,ri— l)j 

fm,n— 1 1 lftm,n) — ^(VVra+l,n VVre,n— l)i 

fm-l,ri{'4 , m-l,n+l ~ 1pm,n) = ^(V'm.n+l ~ 1pm-l,n)- 

Adding these four equations up, one obtains the 5-point scheme (|22|) . 

We are therefore lead to the following change of variables in Zg (see figure 2): 

m — n n + m , 

corresponding to a 7r/4 rotation of the axes and, in the new variables, equation (|22|) is the 
affine self-adjoint five point scheme □ 

Remark. Notice that the potential p satisfies the constrained equation (|12|) 

Remark. Any self-adjoint 5-point scheme © is gauge-equivalent to an affine self-adjoint 
5-point scheme (fTTj), Indeed, if oy^j, is a solution of © with coefficients a^ u , 6^ and c^ )U , 
then 

satisfies the affine equation with coefficients 
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(m-l,n+l) <=> 



l,v) 



(m+l,n+l) <=>((j,,v+l) 



(m-l,n) 0- 




(m+l,n) 



(m,n) <=> (\i,v) 



, / (m,n-l) \ , 

-O— — * 

(m-l,n-l) <=>((j,,v-l) (m+l,n-l) 





<=> (|X+1,V) 



Figure 2: The white (odd) points are eliminated taking a suitable linear combination of 
four adjacent 4-point schemes. What remains is a 7r/4 rotated 5-point scheme on the black 
(even) points. 

2.2 Prom the sub-lattice to the lattice 

The results of the previous section 2.1 leave open the question if all the solutions ip : Z 2 — > C 
of the 5-point scheme (|2*2*|) can be extended to solutions of the 4-point scheme (|2*U)) on the 
whole lattice. The answer is affirmative and the construction is very simple. 

Suppose one knows a solution ip : Z 2 — > C of the 5-point scheme (|22|) for some given 
coefficients a, b : Zg — > R. Define the function / : Z 2 — > R on the whole lattice by the 
inverses of equations (|2*5|) and (|25[). and choose the value of ip at one arbitrary odd point 
of the lattice. Then, using the 4-point scheme (|20[) and the known values of ip at the even 
points, one can construct uniquely a solution ip of the 4-point scheme (|2Uj) whose restriction 
to Z 2 coincides with ip. 

Summarizing the content of the last two sections, we have shown that any solution of 
the 4-point scheme once restricted to Zg, generates a solution of the 5-point scheme 
(|22j) . Viceversa, any solution of the 5-point scheme ()22|) can be obtained restricting to Z 2 
a suitable solution of the 4-point scheme ()2U|) . 

2.3 Discretization of the Schrodinger equation via the Moutard transfor- 



In this section we show that also the method of "discretization via transformations" |19j . 
when applied to the elliptic Moutard (Schrodinger) equation (jHJ), leads to the 5-point scheme 
thanks to the sublattice approach. Consider the general solution ^ of the Schrodinger 
equation © together with its particular solutions 6 1 and 9 2 . One can easily check that the 



mat ion 




(26) 
(27) 
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satisfy again the Schrodinger equations ([8jl. but with the new potentials 



Denote by a a solution of the system 



Bo ) + ( 63 



■ uu. 



then the functions O 1 ^ and ^ 2( - 1 ^ given by 

<r = 0*01(2) = _0l02(l), 



(28) 

(29) 
(30) 

(31) 



satisfy the same equations as, respectively, Vl^ 2 ) and \I/(i). Finally, the function ^I^ 12 ), 
obtained from the superposition formula 



\j/( 12 ) _|_ \j/ 



^(2) _ 
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satisfies the Schrodinger equation (jHJ) with potential 



■ mi 



(32) 



(33) 



and is simultaneously the transform of \I'( 1 ) via 2 (i), and the transform of via 1 ( 2 ). 

Therefore the superposition principle of the Schrodinger equation (JHJ) is described by 
the 4-point scheme ()32|) of the discrete Moutard type, which we know is not a proper 
discretization of jSJ). This seems to be in contradiction with the general rule that the 
superposition principle of a continuous system provides an integrable discretization of it. 
However, as we know, the sub-lattice approach resolves this problem, since the 4-point 
scheme (|32|) reduces, on its subgrid, to the 5-point scheme (see also section l2~4")l . which 
turns out to be the proper discrete analogue of © 24 . 



2.4 Different gauge forms of the 5-point scheme 

For the construction of the Darboux transformations and in other applications it is useful 
to introduce the r-function of the 4-point scheme (|20|) 

Tm+1 '" Tm ' n+1 = fm,n. (34) 
I'm , n T~m+ 1 , n + 1 

The fact presented below, which can be checked by direct calculation, explains the intro- 
duction of such a potential in a broader context. 

Theorem 3. Consider a lattice : 7L N — > V , dvaiV > N > 3, satisfying the following set 
of linear problems 

7 /'m..(i+l)..(j+l)..n ~~ ' l Pm..i..j..n = 9m..i..j.. n {' t Pm..(i+l)..j..n ~ ' l Pm..i..(j+l)..n) > * < 3i (35) 
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then the functions g 1 ^ i ■ n can be parametrized by the potential T m ..i..j.. n 

ij ^"m..(i+l)..j..n^"m..i..(j+l)..n 

9 m ..i..j..n ~ ~ ~ ' z < Ji 

'm..(i+l)..(jf+l)..n 1 m..i..j..n 

and the compatibility condition of the linear system h'J5\) gives the nonlinear system of 
discrete BKP equations 

T m..i..j..k..n T m..(i+l)..(j+l)..(k+l)..n ~ T m..(i+l)..j..k..n T m..i..(j+l)..(k+l)..n+ 

T rn..i..{j+l)..k..n T m..{i+l)..i..(k+l)..ri ~ T rn..i..j..{k+l)..n T m..{i+l)..{j+l)..k..n = 0) i < j < k. 

Let us define the function 

®m,n = ^ 1pm,n, (36) 

then, by direct calculation using equations (|22j) and one can check that the function 
restricted to the even grid Zg, satisfies the equation 

Vh.^aH-i," + ftfti/Vi." + V"+i^*.H-l + V"*A»>"-1 = c^v®^, (37) 

where 

hm,n = j (38) 

7"m— l.n 



and 



C m ,n = I ! 1 ■ I H ! ( ; 1 ; I • (39) 

1~m,n \ 7~m,n— 1 7"m,n+l / Tm,n'Tm— l,n \ T~m,n+1 T~m,n—1 / 

We call equation (|3*7|) the equal-field self-adjoint 5-point scheme. 
If we define the function 

*m,n = - Tm ' H rpm,n, (40) 

T/T m +l,nTm—l,n 

then the function ^ , restricted to the even grid Zg, satisfies the discrete Schrodinger equation 



^M,^ ,T, i ^M-l.^,Tr l ^A 4 ^ ,T, i ^,1/-!. 



+ -P^-v + + = F^^V, (41) 



where 



r m , n = ,p^, (42) 

7~m+l,n 



and 



^m— l,n / Tm+X,n—1 T~m+l,n+l \ T"m+l,n I Tm,— l,n+l T~m—l,n—l \ / ,„\ 

-Tm,n = I 1 I H I 1 I • (4,5 J 

T~m,n \ ^m,n-l 7"m,n+l / T~m,n \ 7"m,n+l Tm,n—1 / 

Notice the following connection formulas between the equal-field and the Schrodinger gauges: 

= ^ \j.,v = ^ i c a»,^ = ^ -^Mi 4 " (44) 
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Remark. The potential r is defined up to multiplication by functions of single variables m 
and n, which implies that also $ and are not unique. 

Finally we mention a direct consequence of the above definitions of the equal-field and 
the Schrodinger gauges. 

Corollary 4. Given a solution 6 : Z 2 — > C of the discrete Moutard equation 1120]) with the 
t -function T mn> then 

a) 6 restricted to the sub-lattice satisfies the affine 5-point scheme ill]) : 

b) the function 

_ T m,n a , ,_v 

Pm,n — "m,nj V*"J 

restricted to the sub-lattice, satisfies the corresponding equal- field 5-point scheme \3 r ]\) ; 

c) the function 



®m,n — _ @m,n — ^ Pm,ni (46) 



Urn r> — — 

restricted to the sub-lattice, satisfies the corresponding Schrodinger 5-point scheme m\ ) 



3 Darboux transformations of the sub-lattice 

An obvious application of the above construction is that, once a sub-lattice of a given 
integrable lattice is identified , one obtains essentially for free some of its integrability 
properties. Here we show, for instance, the construction of the Darboux transformations 
|24j of the 5-point scheme in the affine equal-field (|37|) and Schrodinger (|41|) forms, 

induced by the transformations of the 4-point lattice (|20[). 



3.1 DTs of the affine 5-point scheme 

We first recall relevant material [2Zj on the Darboux-Moutard transformations of the 4-point 
lattice. 

Proposition 5. Given a solution 6 of the discrete J^-point scheme i2(J]) . then any solution 
if) of the compatible linear system of the first order 

— 7, ! (^m+l,n ), (47) 

"m+l,n 

+ VVn), (48) 

"m,n+l 

satisfies the J^-point scheme 12()]) with the transformed potential 

f _ @m+l,n@m,n+l r (AQ\ 
Jm,n — „ „ Jm,nt 

while the corresponding transformation of the r -function takes the simple form 
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Remark. Notice that the function 

&m,n = ~7. (51) 

is a solution the 4-point scheme Q2UJ1 of ij). 

Corollary 6. One can interpret the transformation as a shift in the third dimension of 
the lattice, and the transformation equations \4^-\4^ as linear problems of the form \2b]) 
involving that dimension. In particular, the transformation rule of the r -function 15(J\) and 
the parametrization °f ^ e potential / m>n in equation 12(J\) satisfied by Q m ,n lead to 

equation 

T~m,n'i~m+l,ri-\-l ^'m,n^'m+l,n+l H "7"ra,n4-l ^"m+l,n ^~m+l,n^"m,n+l)) (52) 

of the form of the discrete BKP equation for N = 3. 

Again, the function ip mn , when restricted to the even (or odd) lattice, satisfies the 
affine 5-point scheme However, the transformation equations (|47l) - (|48|) depend on 

values of the transformation potential 6 on the full lattice. Our goal is to constrain the DT 
of the 5-point scheme to the sub-lattice as well. Obviously, given a solution of the affine 
5-point scheme on the sub-lattice, it can be propagated to the full lattice in the spirit of 
section l2~2*l and then used to construct the transformation. We will show, however, that the 
transformations of the 5-point scheme can be done, in a more elegant way. In particular we 
will constrain the transformation equations to the sub-lattice. 

Lemma 7. Let VVn.n and Q m ^ n be solutions of A20\) . and let ip m ,n be the transformed solution 
constructed via J^7| )- J^<S| ). Then the function ip m ,n, defined by 

^m,n = 20m— l^nVVre— l,n; (53) 

with the functions a m ^ n and 6 m>n defined by equations \23\) . satisfies 

1pm+l,n—l VVre.n — ^m— l,n— 1 ifira— l,n+l VVn,n Gm,n^m—l,n+l) i (54) 
VVfi+l,n+l 1pm,n — 0"m—l,n+l{^rri—l,n+l' l Pm,n &m,n^m— l,n+l ) • (55) 

Proof. Substract equation 1)181) from equation (jlTI) . use the four point scheme (|20(l and 
evaluate the result at (m — l,n — 1) to get equation (|54j) . Similarly, add equation (|47|) 
evaluated at (m — l,n) to equation (|i8)) to obtain equation (|55f) . □ 



Corollary 8. The function ip m , n satisfies equation 



Q-m,?i(VVfi+l,n— 1 ^m,n) ^ ^m-l,n+l('/'m-l,n+l ^m,n)~\~ 
b m ,n{lprn+l,n+l ~ 4>rn,n) + &m-l,n-l (VVrt-l,n-l ~ 1pm,n) = 0, 



with coefficients 



*m,n — , a a ' m ' n — „ a a 

Om—l ; n—l"m,n"m—l,n—l "m- l,n-\-X"m,n"m— l,n+l 



(56) 



(57) 
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¥m,n/ 

--o- 



(3~ 



4 



-P- 



■ ¥m+l,n+l 



^Kn,n 



Vm+l,n 



Figure 3: The sub-lattice reduction of the discrete Moutard transformation. 

Proof. The function VVn,n satisfies equation (I22() with coefficients a mjn and 6 m)n obtained 
from / mjra via equations (J23J). Because the function l/9 m ,n satisfies the same equation as 
ipm,n, then 6m t nipm,n satisfies again equation (|22|). but with new coefficients (compare with 
remark 2 after proposition 0). □ 

The following theorem can be obtained just restricting the above results to the sub- 
lattice, but it can be also proven directly. 

Theorem 9. Let tp^u satisfy the self-adjoint affine 5-point scheme Ml}) and let 0^ )V be a 
particular solution of Ml}) . Then the solution tpav of the compatible linear system 



1^H,v+l — 1pn,v = —O'n-l,v(0fj,-l,u1ptJ,,v — 9n,v1pn-l,v)i 

satisfies the self-adjoint affine 5-point scheme I lip with transformed coefficients 



1 



i 1 ^ 



Q/i— fi,u" /l— l,v 



(58) 
(59) 

(60) 



Proof. The compatibility condition of ()58 [) -(|59 [) is the self-adjoint affine 5-point scheme (|1 1|> . 
Equations ()58 [) -(|59 [) and their compatibility condition imply the self-adjoint affine 5-point 
scheme ((TT|) for ip^ v with coefficients given by (|6U|). □ 

Remark. The transformation l|58 |) -l|59 |) is a restriction of the DT of the general self- adjoint 
5-point scheme given in [21] to the class of self-adjoint affine 5-point schemes. 

Remark. Equations (|58 |) -(|59 |) are obtained restricting equations (|54 j) -(|55 |) to the sub-lattice, 
while equation l|oT)|) is the corresponding restriction of equations (fSTI) . 

It seems that the transformation acts from the even (black) lattice Zg to the odd (white) 
lattice (or vice versa). However, the proper point of view is to consider (see corollary 
EJ) the transformation as the shift into an additional dimension of the lattice, and then the 
transformation acts between black (or between white) points of the lattice (see figure 3). 
There are three other equivalent choices of the transformation on the small lattice which 
give also restrictions of the DT to the sub-lattice. The resulting transformation formulas 
differ from (|58j l -(|59j l in some shifts. The present form has been chosen to be in agreement 
with the earlier formulas of |24| . 
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3.2 DTs of the self-adjoint equal field 5-point scheme 

Let us rewrite equations (|54 |l -(|55 |l using, instead of the fields ip m>n , am,m ^m,n and 6 mjn , the 
fields 3> m ,n> /i m ,n and p m ,n defined by equations (|3l)|) . (|38f) and The natural counterpart 
of the function given by equation (|53")) is the function 

&m,n = i$ m -l,n, (61) 

where Q m ,n is given by the transformed version of equation (|36j) . i.e., 

$m,n = ~ Vm,n = 7 • ( 62 J 

The counterpart of the DT of the affine self-adjoint 5-point scheme described in theorem 
will be the following theorem on the DT of the self-adjoint equal field 5-point scheme. 

Theorem 10. Let and p^^ be solutions of the self-adjoint equal field 5-point scheme 
\'JT\j , then the solution of the following linear system of the first order 

' ( 63 ) 

Pw+ih^p + i<&^ v+ i - p^uh^u^^u = -h^ u (p^i-i^^u - P^v^^~i,v) , (64) 
satisfies the self-adjoint equal field 5-point scheme with coefficients 

K = (65) 

P)J,-l,V-l 

and 

c^ v = hwPw + + — s — — + — . (66) 

\Pn,u-l Pn-l,vJ Pn-l,v-l \"-n-l, v Pn-\ 

Proof. The theorem can be verified by direct calculation. To do its derivation within the 
sub-lattice theory (the form of the transformation equations has been explained before the 
formulation of the theorem), notice that the function hav is restriction of 

h m ,n = hm—l,m with h m ,n = ~ — ■ (67) 

T~m—l,n 

Also the transformed coefficient c„ v is the restriction of 

Cm,n — Cm— l,n.) (68) 

with Cm,n given by equation (|39|) . with T m n instead of T m n . In putting it on the sub-lattice 
one uses equation (j2U|) satisfied by 6 m ,n- d 



Remark. In terms of the function 



Pw = 7 — , (69) 

Pn,i/ilfi,u 
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formula Iffifijl takes the form 

+ ^n,vP^-\,v + hfj_ iU+ ip^ )U+ i + h^vp^v-i j , (70) 

Pfi,v 

i.e., the function pn U is a solution of the equation satisfied by & m ,n- This is the 5-point 
analogue of the fact that the function Q m ^ n = V^m.n satisfies the transformed equation (|20[). 
Indeed, the function is the restriction 

Prn,n — Pm—X,n C^^-) 

of the function 

~ _ r m,n g 

Pm,n — ~ "m,n: v' / 

obtained from mn in the same way as $ mjn is obtained from t/j mn . 
3.3 DTs of the Schrodinger 5-point scheme 

In what follows we will obtain the DT of the 5-point Schrodinger scheme |24j from the sub- 
lattice theory. In order not to repeat the full procedure of the previous two subsections, we 
will use the connection formulas (|44[) between the equal-field and Schrodinger gauges and 
the form (|46|) of the particular solution Q^ yU of the Schrodinger scheme. The transformed 
potential T^u can be obtained from the connection formulas (|44|) and the potential ha,v of 
the equal- field scheme given by equation 1)65 Jl . i.e., 



rjU - ~ — - r^^r^-x^-i " M n 1 '" 1 ■ (73) 



The transformed solution ^^u, because of (|44|1. is related to by 

f 

The transformed potential Fa v can also be obtained from the connection formulas 1)44(1 



1 l'-'/i,!/— 1 L /i—l,u y ^'ii—l,u V^A 4 — 1,1' ^z^,^ — 1 , 

(75) 

Rewriting theorem llUl in the new fields, we obtain the DT of the Schrodinger 5-point scheme. 

Theorem 11. Let \P ^ v and 0^ be solutions of the Schrodinger 5-point scheme 141\) ; then 
the solution ^ ^ v of the following linear system of the first order 

^±^6^!,^+!,, - ^e^*^ = (e^-i^w - e ft „Vi) . ( 76 ) 
E^e^+^+i _ E^e^*^ = -Epi^ (e rl) „v - e M ,^ M -i,„ ) , (77) 

where T^u is given by equation {7ty , satisfies the 5-point Schrodinger scheme with the 
coefficients r M)J/ and the coefficients given by \75\j . 
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Remark. The analogue of the self-transformed solution Puv, given by equation (|fi9|) . of the 
equal-field scheme is the solution 

of the transformed Schrodinger scheme. In terms of the formula (|75jl for the trans- 
formed potential F^ v takes the natural form 

K» = a~ f J^e/H-i,, + %^<W + r^QwM-i + ^Q^-il • (79) 
4 Superposition of the DTs 

Using the superposition principle of the DTs of the discrete Moutard equation [27], we will 
obtain the corresponding superposition principles for the DTs of the self-adjoint 5-point 
schemes in the equal field gauge and in the Schrodinger gauge. 

Proposition 12. Given solutions '4 , m,n, ^m,n ana " ®m,n °f the J^-point scheme \2U}) . denote 

byipm]n andipm]n the transforms respectively ofi\) m ^ n via9m n and6^ nn . If o~ m ,n is a solution 
of the compatible linear system 

0~m+l,n ~ 0~m t n = 9 m ,n®rn+l,n ~ @m,n@m+l,n' (80) 
0~m,n+l ~ 0~m,n = @m,n@m,n+l ~ @m,n@m,n+l> (^-Q 



then 6m\n and 9m n , given by 



m# =_fl2(i) fl i = (82) 



'm,n m,n " m,n m.n 



are solutions of the ^-points schemes satisfied respectively by tpm,n and^m]n- Moreover, the 



function iplufn given by 



nl q2 

V'min — ^rn,n = 1 ~ (VmiL ~~ V'm^n)' (83) 



a 



in.it 



is simultaneously the transform of ipmin via Qm,n and the transform of ipmln via 9m,n ■ The 
corresponding transformation of the T-function is given by 

T m,n = CF ra,n'Tm,n- (84) 

Our goal will be to derive the corresponding Bianchi superposition principle for the 
5-point equal-field scheme, and then for the Schrodinger scheme. First, we will rewrite 
equations l|8U |) -(|81j l in terms of the transformation data 

(4n,n = ~ ' ^m,ni j = 1 5 2, (85) 

' m+l,n 
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of the equal-field gauge (see equation (j45j0 . and we will put them on the sub-lattice. 
From equations (|8U j) -(|81 j) and the discrete Moutard equation (|2U|) we have 

_ . Tm,nT~m—l,n—l ( a \ q2 _ Q 2 nl \ /of?\ 

C"m,n-1 l,n — 1 \"m,n"m— l,n— 1 "m,n"rn—l,n—l) ; l 0D J 
im—l,n'm,n—l 

_ _ _ ^ r m,raTm-l,ra+l fl 2 /i2 zil \ /o 7 ^ 

Vm,n+1 »m-l,n — L ~ ~ (, t7 m,n t7 m-l,n+l a m,n a m-l,n+l) • I ' J 
'm— l,n 'm,n+l 

If we introduce the function 

S m ,n = i<?m-l,n (88) 
then the above equations can be rewriten in the form 

S m +l jn -l — S m)n = /irn,?i {Pm : nPm-l,n-l ~ Pm,nPm-l,n-l) > (89) 
S m +l 5 n+l ~~ 5]»n,n = ~h m ,n {Pm,nPm-l,n+l ~ Pm,nPm-l,n+l) • (90) 

The functions <3?m^n, constructed according to equations l)61 |) -l)62 [) 

Q3 ,/.(?') 

*%n = * ^l" 1 - 1 ^ , J = 1,2, (91) 

are the transforms of $/i,^ via pjj,,u, when restricted to the sub-lattice. Consequently, the 

j(k) 

transformed data p^ u , j 7^ k are restrictions of the functions 



/^S = « ' ,m " fc 1, \ m " 1,W . ^ (W) 



TO— 1,71 ) 

which, due to equation (|82|) . are related to E m)Jl by 

z, J2 1(2) _ _l 1 „2(1) _ s .qo^ 
""m,nl J rn,nl J m,n ~ ""rn,nl J m ,nt J m,n ~ ^-"m,n- \ ua J 

Using equations (|47| ) -(|48 |) we modify the superposition formula ()83|) in order to include 
functions Vw-l,n an ^m-l,n insteacl of ipm,n and V>m?n 

C"m,nV'm,n — °77i— l,7i' i / , 77l,7i = @m,n@m—l,n l Prn—l,n ~ @m,n@m—l,n' l l J m—l,n- (94) 

Finally, applying twice equations (|61j) - (|62j) . we obtain that the function 

.(12) 
m-2. 



$(12) _ T m-2,n ,(12) , q _x 



(12) 

gives the desired superosition , when restricted to the sub-lattice. 

Putting all the above considerations together, we obtain the Bianchi-type permutability 
theorem, which can also be verified by direct calculation. 
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Theorem 13. Given solutions §^ u , Pjj, v and p 2 „ of the 5-point equal field scheme j37\) , 

denote by $^ and the transforms respectively of via and pi,v tf^w is a 
solution of the compatible linear system 

^ft+l,u — = h^uip^uP^v-i ~ p\,vP{i,v-\)i (96) 

Sju.iz+i - S/i.i/ = -h^vip^vp^-i^u - p\,vPii-\,v)) (97) 

i/ien i/ie functions p]ffi and p'fffl given by 

h o^o 2 = -h o 2{l) o l = £ (981 

are solutions of the 5-points schemes satisfied by $J^l a^rf correspondingly. Moreover 

(12) 

i/ie function <&)i,v , given by 

, 12 

(12) _ h-wP^Pw . (2) (1) , , , 

Mi" 

is simultaneously the transform of 3*^1 via transform of <l>j^|t uia /O^^ . 

Corollary 14. Equations 16 5\) and \6b)) . together with the above theorem, imply the trans- 
formation law of the coefficients of the doubly transformed equal field scheme 

, (12) _ S^+i^+i 



„( 12 ) 

3 /i+l,v+l 



r 2 

f^+l ; H-l 
S/i+l^^Ai^+l 



'A*i^ 



/l 2 

ll^u / 1 2 1 2 \ 

+ ^ -L ~ \P n,u-lP n~l,u ~ Pfi-l,vPn,v-l) + 



+ 



2 _ 1 2 \ 

P)jL,v+lPn+l,v Pfi+l,vPfj,,v+lj 



■ (ioi) 



Finally, let us formulate the corresponding superposition principle for the Schrodinger 
scheme. It follows just from the connection formulas Q44D and (|46j) applied to theorem 1131 

Theorem 15. Let ^,1/? anc ^ ^ e solutions of the 5-point Schrodinger scheme \41\) , 



denote by and the transforms respectively of^^^ via 9* and O 2 ^, and let F 
denote the first coefficients of the corresponding equations, i.e., 



(i) 

At, i/) 



r (i) 



A' — 1,^—1 -p -p „• _ i o 



is a solution of the compatible linear system 



(102) 



r,. 



)J,,L> 



^M^+l — ^A*,^ ~~ p ~(®li,v®fi-l,v ®aU®M~V )' 
1 A^ 



(103) 
(104) 



17 



then and &f$ , given by 



r (2) p(l) 
B 1(2) B 2 _ _ 2(i) i LE£. - V 

1 Mi^ 1 A*." 



(105) 



are solutions of the 5-points schemes satisfied respectively by Vpj^j, and W^jv- The function 

(12) 

^fi,u , given by 



VhV+l r (12) ^(12) p ^ _ e ^ e ^ fr (2) ,t,(2) r (l),T,(l) 



(106) 



where 



,(12) 



S^t+1,!/+1 

«2(1) 



(107) 



is simultaneously the transform of ^f$ v via Qf^J) and the transform of via 

(12) 

satisfies the 5-point Schrddinger scheme with the coefficient T^J given in \lUT\j , and 



with the coefficient F^ u given by 



F, 



(12) 



,(12) 



+ 



V p2 



r 2 



s^+i^+ir^+i^r^^+i 



(0^+i0"+i,, - 0i+i,,0^+i) ) • (108) 



5 Algebro - geometric solutions 

In this sections we construct algebro - geometric solutions of the 4-point scheme 1)20(1 and, 
due to the results of the previous section, of the 5-point scheme (|22 [l .(|ll |l . 

The direct and inverse periodic transform for a generic 4-point scheme (1121) was devel- 
oped in ^Hj. To obtain algebro-geometric solutions of a specific 4-point scheme we have 
to impose suitable constraints on the spectral data. As we shall see, the constraints which 
give rise to the 4-point scheme (|2"U|) are analogous to those introduced in [3U] in the study of 
the 2D continuous Schrodinger operator. Constraints of such type were first introduced in 
[S| in the theory of reductions of 1+1 systems. In the theory of discrete systems, analogous 
reductions were used in [Q to characterize orthogonal and Egoroff quadrilateral lattices. 

We remark that algebro-geometric 5-point schemes were also studied in ^S]. The main 
difference between our 5-point schemes and the schemes studied in ^Hj is the following. In 
our case, the full Bloch variety for zero energy is constructed and the eigenfunctions with 
other energies are not incorporeted in the spectral transform. Moreover, only self-adjoint 
schemes are constructed. In |18j . generic schemes with the following property were studied: 
for each eigenvalue, a finite number of eigenfunctions are explicitly constructed. 

The standard finite-gap construction JB] for a generic 4-point scheme is based on the 
following spectral data. Assume that we have: 

1. A compact, regular, connected Riemann surface T of genus g. 
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2. / + 1 points Ri, Ri+\ in T - the normalization points for the wave function 
(denoted by "if). 

3. / + g points 71, . . . , Ji+ g in T - the divisor of poles of the wave function. 

4. A collection of points Pf , . . . , P^, P-f, . . . , P^, Qf, . . . , Q~j^, . . . , Q~^, where 
M, N are arbitrary positive integers. 

From the Riemann-Roch theorem it follows that, for generic data, there exists an unique 
function 'I' (7, 777, 77), where 7 G I, 777,77 G Z , 1 < m < M, 1 < n < N, with the following 
properties: 

1. 'I' (7, 777, 77) is a meromorphic function of 7 in T. 

2. ^(7,m,n) has at most first-order poles at the points 7^, k = 1, . . . , g + /, P£ , k = 
1, . . . , m, , k = 1, . . . , n and no other singularities. 

3. <I>(7,m,ra) has at least first-order zeroes at the points PZ , A; = l,...,m, A; = 
1, . . . , n. 

4. fy(R k ,m, n) = 1, k = 1, . . . , I + 1. 

The discrete time shift m — > m + 1 corresponds to adding one extra pole P^+i an d one 
extra zero to the wave function. 

Let us check that function ^f(^f,m,n) satisfies the 4-point equation 

^(7, m+1, n + l)+«i(m, 77) ^ (7, m + 1, 77) +02(777, 77) ^(7, m, 77+ 1) + 03 (777, 77) 'I' (7, m, n) = 0, 

(109) 

in which the 7-independent coefficients ai(m, n), a2(m,n), 03(777,77) are defined by the 
following formulas: 

01(777,77) = — hm — — — (HO) 



m + l 



^(7, 777 + 1, 77) 



a 2 (m,n) = - hm — — — — (111 

0:3(771, n) = —1 — 01(777, rt) — 02(771,, 77). (H2) 
Indeed the left-hand side of (|1U9|) has the following properties: 

1. It is a meromorphic function of 7 in T. 

2. It has at most first-order poles at the points 7^, k = 1, . . . ,g + I, P^, k = 1, . . . , m, 
Q~£, k = 1, ... ,77 and no other singularities. (Conditions Q11U|) . (jlllj) mean exactly 
that the poles at the points P^ + i, Qn+i vanish). 

3. It has at least first-order zeroes at the points P k ~ , k = 1, . . . , m, , k = 1, . . . , n and 
R k , fc = l,...,Z + l. 

From the Riemann-Roch theorem it follows that, for generic data, this function is identically 
equal to 0. 
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Remark. In the paper JH] the points P£, PZ : Qf, were generic. The special case in 
which all the P^'s coincide: P£ = P + for k = 1, . . . , M, as well as all the PjT's , Q^'s, 
Qh's: P^ = P~ , k = 1, . . . , M, = Q + , = Q~, k = 1, . . . , N was also discussed in 
|18j . In this case, the discrete time shift m — * m + 1 corresponds to increasing the order 
of the pole at the point P + and the order of the zero at the point P~ of one, and the 
Abel transform in the ^-functional formulas (see below) becomes constant. This last choice 
results in an essential effectivization of the explicit formulas, but it corresponds to a more 
restricted class of potentials. For example, if the potentials are periodic in m with period 
M, the Floquet multiplier at the point Qn+i can De calculated by the following formula 

«i(n) = (-l) M TT a3(m,n l 
- LA 01(771,77) 

m=l 

If all points Q~ coincide, then «i(n) should not depend on n. 

Since we do not like to impose such a restriction on the class of our potentials, we assume 
a generic collections of poles and zeroes. 

5.1 Discrete Moutard reductions and constraints on the spectral data 

Let us describe the reductions corresponding to the affine 4-point scheme ()20|) (may be with 
a complex potential f(m,n)). 

Lemma 16. Assume that T possess a holomorphic involution a with exactly 2 fixed points 
R + = R\, R_. Assume that the spectral data have the following symmetry with respect to 
a: 

1. Exactly one normalization point R\ = i?+ is used (I = 0). 

2. There exists a meromorphic differential 0, with 2 first- order poles at the fixed points 
R + , R- and 2g zeroes at the points 71, . . . , j g , aji, . . . , aj g . 

3. aP+ = aQ+ = Q k for all k. 
Then 

0:3(771, n) + 1 = 0, 01(771, 71) + 02(771, 71) = (113) 
and, consequently, $(7,771,71) satisfies the J^-point scheme \2U\) with 

fl \ ■ i \ -v ^(l,m + l,n + l) -v $(7,m + l,n + l) , 11/n 
r(m, n) = 101(771,71) = —1 hm ; r — = 1 lim ; r — . (114 

K ' U ' 7-P+ +1 *(7, m + 1, 71) 7 ^Q+ +i $(7,771,71 + 1) V ; 

Proof. By analogy with consider the following form: 

0$(7, 771, 7l)$(o"7, 771, 7l). 

This form has 2 first-order poles at the points -R+, R- and no other singularities. Therefore 
res(0$(7, 771, ?i)$(o"7, 771, n), 7 = -R+) = — res(0$(7, 771, ti)$(o"7, m, 71), 7 = R-) 
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Taking into account that 

res(f2, 7 = R + ) = — res(fi, 7 = R-), 

we obtain 

ty 2 (R_,m,n) = fy 2 (R + ,m,n) = 1. 

Let us show that 

*(i?_,m,n) = (-l) m+n . (115) 

If m = n = we know that ^(7,0,0) = 1 (from the uniqueness argument and taking into 
account that the constant 1 satisfies all the analyticity constraints). Therefore ^f(R-, 0, 0) = 
1. Consider the function ^(i?-, m + l,n) as a function of the point P^ +1 (we use here that 
^m+l = ^m+i) 1 By construction this function is meromorphic in P^+i- Since the surface 
T is connected, then ^(i?_,m + l,n) does not depend on P^+i- Assume that -P^ +1 be 
located very close to i?_ and denote by z a local coordinate in a neighbourhood of 
z(R_) = 0, such that ctz = — z. Then 

2 + P + 

V(z, m + 1, n) = _ m, n) + o(P+ +1 ); 

therefore 

V(R-,m + 1, n) = -$>(R-,m,n). 

Substituting (|115|) into Q1U9|) and taking into account l|112j) . we obtain (|113|) . This 
completes the proof. □ 

Having in mind applications in discrete geometry and in discrete complex functions 
theory, let us formulate sufficient conditions which guaranty that the potential f(m,n) be 
real. 

Lemma 17. Assume that all the constraints of Lemma \lb\ on the spectral data be fulfilled. 
Let the spectral data satisfy the following additional restrictions. There exists an antiholo- 
morphic involution r onT with the following properties: 

1. t commutes with a. 

2. tR + = 

3. The points P^ , P^ , , are fixed point of r. 

4- The set of divisor points yj~ is invariant under r (but r may map one divisor point 
into another). 

Then the potential f(m, n) and the function ^(j, m, n) have the following reality properties: 



f(m, n) = f(m, n), *(t7, m, n) = (-l) m+n ^( 7 , m, n). (116) 

In particular, if 7 is a real point ofT: tj = 7, then ^(7, m, n) is either real or pure 
imaginary: 

\l/(7, m, n) £ M. if m + n is even , , 

^(7, m, n) £ itik. if m + n is odd. 
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To prove the Lemma, it is sufficient to point out that the function (— l) m+n ^(r7, m, n) 
satisfies all the analytic constraints characterizing \&(7, ret, n). Therefore these two functions 
coincide. The reality of f(m,n) then follows immediately from (|114|) . 

Remark. The reality constraint on the eigenfunction ^(7, m, n) naturally distinguishes be- 
tween even and odd sub-lattices. 

In Section 2 it was shown on the algebraic level that any solution of the 4-point scheme 
(|2fl|) generates a solution of the 5-point scheme (flip. In this section we give an analytic 
proof of this statement. 

Consider the following meromorphic 1-form 

0(7, m, re, rh, re) = 0(7)^ / (o"7, m, re)\&(7, rh, re). 

Assume that the differential £1 used before to define the constraints on the spectral data 
has the following normalization: 

= -, res(f},7 = R-) = --. 

Let us calculate the residues at the poles of Cl{p/,m,n,rh,n) for \m — m\ < 1, |re — re| < 1. 
It is convenient to write the answer in graphic form. We denote our Riemann surface by an 
oval, and we assume that our distinguished points are located at the following positions: 



res(ri, 7 = 




In the next picture we show the poles of ^(7, ret, re, rh, h). We use the following notations: 
+ denotes a pole with residue — denotes a pole with residue — |, x denotes a pole with 
a residue different from ±5, o denotes a zero. The parameter n is written in the left column, 
the parameter rh is written in the bottom row. The residues are written next to the symbols 
x. 



-if'Cm-Ln) 



n+1 



n 



ir(m-l,n) 



-1 






-if(m,n) if(m,n) 



n-1 



-if(m-l,n 




if(m-l,n-l) 




-1 




irtm,n-l) 



if (m,n-l) 



m-1 



m 



m+1 



22 



Let us explain how the residues are calculated. First of all, for to = m, h = n we 
have only 2 poles and the residues are opposite, therefore we have an identity (which was 
used when we proved the Lemma llfljl . Ifm = m±l,n = norm = m,n = n ± 1 we 
have 3 poles. At R + and R- the residues are +|, therefore the 3rd residue is —1. If 
\rh — to | = |n — n\ = 1, the residue at R + is equal to the residue at R- is equal to — ^, 
therefore we have 2 more poles with opposite residues. For example, let us calculate the 
residues for to = to + l,n = n + 1. From the 4-point scheme (|20|) it follows that 

res(0(7, to, n,m + l,n + 1),7 = -P+ +1 ) = if(m,n) res(f2(7, to, n, to + 1, n),7 = 

therefore 

res(fi(7, to, n, to + 1, n + 1), 7 = = —if(m, n). (H8) 

All the other residues can be calculated exactly in the same way. 
To prove that the function ^(7, to, n) satisfies the 5-point scheme 

7 ^ r )(^( 7 , to + 1, n + 1) - *( 7 , m, n)) + 

f {m -\,n-i) (^(7. to - 1, n - 1) - ^(7, to, n))+ (ng) 
/(to, n — l)(^(7,TO + l,n — 1)— ^(7, m, n)) + 
f(m — l,n)(^ r (7, m — l,n + 1) — ^(7, m, n)) = 

it is enough to verify that the left-hand side of ()119|) has the poles and zeroes prescribed 
for \P(7,to, n) and, moreover it is equal to zero at the points R + , R-. Therefore, by the 
Riemann-Roch theorem, it is zero. 



5.2 Algebro-geometrical solutions 

Riemann surfaces with the constraints described above can be constructed in the following 
way. We start from the "vacuum" solution, corresponding to the Riemann sphere g = 0: 

If all Pfc's coincide as well as all Qfc's: Pk = P, Qk = Q, the above vacuum solution 
reduces to the so-called "discrete exponential" 

discussed in [1]; if Q = iP, then /(to, n) = 1 and (|12U|) solves the Cauchy- Riemann equation 
Q. 

The involutions a and r are given by: 

a\ = —A, tA = ■=■ 
A 

and the Riemann sphere with all its distinguished points is drawn here: 
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By analogy with the continuous case, discrete holomorphic polynomials can be con- 
structed from the discrete exponential (|12U|) in the following way (if Q = iP): 



p^(m, n) 



1 d?e{m, n; l/£) 



The first 3 examples read: 



,(0) 



1, p 



(i) 



m + in, p^ = (m + in) z 



d( 3 



,(3) 



(121) 



(m + in) 3 + (m- in) /2. (122) 



We observe that the discrete holomorphic polynomials coincide with the continuous mono- 
mials z 3 (assuming z = m + in) up to the order 2. 

A non-trivial Riemann surface can be constructed by attaching handles to the above 
Riemann sphere keeping the above symmetries. Two examples are shown below. 



real ovals 





The figure on the left illustrates the so-called M-curve (the number of real ovals, i.e. the 
ovals formed by the fixed points of the antiholomorphic involution r, is equal g + 1. g + 1 
is the greatest possible number of real ovals. In this particular case g = 4). The figure on 
the right illustrates a curve with g = 4 and only one real oval. 

Remark. In the theory of the 2-D continuous Schrodinger operator the M-curves play a 
distinguished role. In particular the operators generated by the M-curves are non-singular 
and strictly positive [SO], [Tl]. A complete classification of finite-gap data corresponding to 
nonsingular solutions was obtained in |23j . The study of analogous properties of the 4-point 
and 5-point schemes will be the subject of future investigation. 
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It is well-known that function ^(7, m, n) can be written in terms of Riemann ^-functions 
in the following way: 



VP (7, m, n) = 



9 A( 7 )+ E E (A(Q»)-^(0+))- £ A^)-* 



fc=i 



fc=i 



fc=i 



B 



9[A{ 1 )-Y J A{ lk )-K 
k=i 



D 



(123) 



9^A(R + )-^l( lk )-K B^j 



B 



9 A(R + ) + E (A(p-)-A(P+))+ E (A(Q^)-A(Q+))- EMlk)-K 
\ k=i k=i k=i 

(ml. n 1. 

E J n ^ p k' p k) + Y, J n ^Qt,Qk 
k=l R+ k=l R+ 

Here we have used the following notations. a k , b k are the basic cycles in T, a k ° h = 8 k u 
a k o ai = b k o bi = 0. uji are the basic holomorphic differentials such that 



&i — Sij, 



(124) 



$7(7, P, Q) are meromorphic differentials of the third kind with 2 first-order poles in P, Q, 
with residues —1 and +1 respectively and zero a-periods: 

res(n(7,P,Q), 7 = P) =-1, res(n( 7 , P,Q), 7 = Q) = 1, / fi( 7 ,P,Q) = 0. (125) 

^(7) = (^1(7), ..,.A ff ( 7 )) denotes the Abel transform: 

7 

A fc ( 7 ) = y w fc , k = l,...,g, (126) 

where the starting point V of the Abel transform can be chosen arbitrarely. The Riemann 
^-functions are defined by the following Fourier series: 



9{z\B)= E ex P< 

mi,...,m g 

where B k i are the 6-periods of the holomorphic differentials: 



7ri B k jm k rrij + 2ttz z k m k 

kj k 



(127) 



oji = B 



3* 



(128) 
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(in the definition of the ^-function we use the same normalizations as in j2Hj ; other classical 
books, like [TSJ, use the different normalization corresponding to a-periods equal to 2m). 

Due to the symmetry of the spectral curve imposed by a, the Riemann ^-functions can 
be expressed as bilinear combinations of Riemann and Prym ^-functions of genus g/2 \12\ . 
This is analogous to the continuous Schrodinger case |30] . 

The simplest non-trivial example corresponds to an M-curve with 2 handles. 




In this case, due to the symmetry imposed by a, the answer can be expressed in terms of 
genus 1 ^-functions. It is convenient to choose a basis of cycles with the following symmetry 
(see [12]): 

aa\ = —a 2 , ob\ = —b 2 . (129) 
Denote by u and Cj the holomorphic differentials such that 



-■j.:. J) u) = J> a — I 

ai ai 



and let 



Then 



and 



rj = <p u>, ;/ :/) ^. 

bi bi 



U! + 00 
^1 = — , 0J 2 



<jJ — UJ 



B 



11 



fj + rj 

B 22 = — ^ — , B12 = B21 



7] — 7] 



d( Zl , z 2 \B) = 6(zi + z 2 \2r ] )9{z 1 - z 2 \2fi)+ 



9(zi + z 2 + rj\2i])9(z 1 - z 2 + t)\2t)) exp 



m ( + 2zi 



(130) 



(131) 



(132) 



(133) 



Let us consider the following explicit example in which the Riemann surface V is defined 
by the equation: 



fJ(A-e fc ), e 7 _ 



-efe, 



(134) 



k=l 



R+, are the infinite points /i ~ A 3 and 11 ~ —A 3 respectively, with A ~ 00, a : (X,fi) 
(—A, — fi), t : (A, /u) — > (A, — /x), u and lu are defined by: 



, d\ „ , XdX 
u = d\ — , to = d 2 . 



(135) 



where the parameters d%, d 2 are fixed by the normalization condition (|13U|) . The a and b 
cycles are shown in the figure below, the cuts are along the intervals [e\, e 2 ], [e^, e^], [e^, e^], 
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the solid lines correspond to the sheet containing R + , the dashed lines correspond to the 
sheet containing 




The divisor is located on the ovals a±, 02, the distinguished points P£, PZ ', Qt, Q^T , 
are located on the oval lying over the interval [e% , . 

If the starting point of the Abel transform is chosen at one of the branch points of a 
hyperelliptic Riemann surface, the vector of Riemann constants can be written in a very 
simple form (see ^U). If = e 6> then 



K 



1 1 B12 

2 < 2 

1 1 £22 

2 2 



(136) 



5.3 Some applications 

As it was already mentioned in the introduction, the 4-point and 5-point schemes Q2U|1 and 
(jlljl are relevant in the theory of discrete complex functions, as well as in the theory of 
integrable discrete geometries. 

Solutions ^(7, m, n) of equation (|2()jl provide examples of discrete holomorphic functions 
and, through the change of variables (|25|). generate discrete harmonic functions on the sub- 
lattice. Let us concentrate here our attention, for instance, on the discrete analogues of a 
polynomial of degree 2. 

If g = 2, the corresponding discrete holomorphic function is constructed from the 
^(7, m, n) function (|123|) . (|127j) associated with the hyperelliptic curve Q134j) . via the for- 
mula: 

lim ([y(i,m,n) + y(aj,m,n)-2ty(R + ,m,ri)]—). (137) 

Due to the positions of the distinguished points P + , Q + ', the corresponding 5 = 
discrete holomorphic function has the following form: 



sft(VlV 2) (™,™)) =&(e i **(m + in) 2 ) . 



(138) 



The graphs of the real parts of the discrete holomorphic functions (|138|) and (|137|) are 
shown in the figure below. 
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y 



y 



g = 



For g = 2 the spectral data are: e t = -3, e 2 = -2.02, e 3 = -2, 71 = (2.6,^16.56635904), 
72 = (-2.6,^16.56635904), P + = (-l,iV73.9296), Q + = (—v/3, -i\/6.4824). 

As it was shown in [25], solutions ^(7, fM, v) of equation (jllj) allow to construct quadri- 
lateral surfaces in IR 3 , i.e., discrete surfaces whose elementary quadrilaterals are planar, in 
the following way. 



Consider the vector N, 



L 2 — ► IR 3 whose components are three independent real 
solutions of the self-adjoint 5-point scheme Then Na V is the normal vector of a 

quadrilater surface : Z 2 — > R 3 defined by the following generalization of the Lelieuvre 
formulas: 



^(J,ffJ,,U 



(139) 



Here we consider examples of quadrilateral surfaces constructed using the discrete ana- 
logues of complex polynomials of degree 1. If g = 0, we consider the following normal 
vector: 

N$l = [25ft (jfrfa -u + i(ji + 1/))) , 29 (e^Ou -u + i{fi + , 1 
Analogously, if g = 2, the normal vector is defined by 

45= 



(140) 



lim ([*( 7 ,^-i/ jM + i/)-1]A) ,9f lim ([*( 7 , p - v, n + v) - 1] A) , 1 , 

(141) 
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where ^ is the wave function (jl23j) . (|127j) associated with the hyperelliptic curve (| 134)1 . 
Then the two quadrilateral lattices constructed, via the embedding ()139j) . by the normal 

vectors Nf^l and N^l, are shown in the figures below. 




9=0 9=2 
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